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Abstract 



We consider heavy quark and antiquark off-shell spinorial currents with emission of 
an arbitrary number of gluons of positive helicity. From this results we calculate the 
corresponding on-shell scattering amplitude and the quark-antiquark vector current. Then, 
we show that in the heavy top quark effective theory the holomorfic component of the Higgs 
— > qq + n-gluon amplitude vanish for helicity conserving configurations. 
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1 Introduction 



Production of jets at large transverse momentum with respect to the beam is the typical signature at high 
energy hadron colliders, both for signal and background. The LHC will produce an enormous amount of 
heavy quarks. Indeed it can be considered a heavy quark factory. Since heavy quarks will appear associ- 
ated to jets, the calculation of multipartonic scattering amplitudes with heavy fermions is mandatory in 
order to provide accurate phenomenological predictions. But the complexity of perturbative calculations 
increases exponentially with the number of jets, in such a way that the usual expansion in Feynman 
diagrams becomes prohibitive. Recursion relations within the helicity amplitude formalism ID El 13 
have been proven to be an elegant and efficient tool to calculate multipartonic scattering amplitudes, 
and might allow to overcome this problem. Together with the calculation of the scattering amplitudes, 
more efficient methods to integrate over the intricate multipartonic phase- space are also required for the 
development of Monte Carlo event generators. 

Recursion relations have been extensively used in the literature at tree [4, 5 1 and one-loop level 
in the past. Witten's idea of a duality between supersymmetric Yang-Mills and topological string theories 
in twistor space JSJ has reopened the interest in this issue. Since then, a new method for the evaluation 
of scattering amplitudes in gauge theories has been proposed [9|. It is based on the recursive use of off- 
shell Maximal Helicity Violating amplitudes (MHV) [ 10 1 . Recent works have accomplished interesting 
progress since the original formulation, and the method has been refined by introducing more efficient 
recursion relations lfTTTfT2l . and extending this approach to the one-loop level lfT3lfT4l . Although ini- 
tially formulated for massless particles the MHV rules have been generalized to include heavy particles. 
Within this formalism recursion relations for massive scalar particles at tree-level have been introduced 
in Ref. ifTol and for vector boson and fermions in Ref. [ 15 1. Those recursion relations have been used by 
the authors of Ref. ifTTI to calculate multigluonic amplitudes with heavy scalars. 

Due to the phenomenological impact of heavy quark production at LHC, we analyze in this paper 
the potentiality of the Berends-Giele recursion relations @ to evaluate several helicity amplitudes with 
heavy quarks and an arbitrary number of gluons. We obtain compact expressions for some helicity 
configurations. We find that this method provides efficiently valuable information for the calculation of 
scattering amplitudes with heavy particles, in a complementary way to the MHV approach. 



2 Scattering amplitudes with one quark-antiquark pair 

Spinors for massive fermions can be constructed from two null vectors |[THlfT^ll2l^ . The decomposition 
is however not unique, and the freedom to choose such vectors makes the helicity amplitude method 
very powerful also for massive particles. We shall consider processes where quark-antiquark pairs are 
produced. We denote by p± and p%, with p\ = p\ = m 2 , the end point four- momentum of the fermion 
line. In terms of two light-like vectors {p\= p\ = 0) the outgoing fermion-antifermion four-momenta 
can be written as 

Pi = -l^—Pl+^—P2i 

m l-(3 H l + ^ 

P2 = —FT- K + P2. (!) 
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where (3 = y 1 — Am? / si 2 is the velocity of the quark, and s i2 = (pi +P2) 2 - This momentum transforma- 
tion is of course not unique. The advantage of our choice of normalization is that it preserves momentum 
conservation since trivially pi + p 2 — Pi + P2, and hence s 12 = (pi + P2) 2 = s 12 . Furthermore, in the 
massless limit we have: p x — > p\ andp 2 — * Pi- The inverse transformation reads: 

~fi _ 1 + P fx 1-P fj, 

P2 ~ ^ Pl+ w P2- (2) 

Then, we use p\ and p 2 as reference null vectors to define the quark- antiquark spinors: 

u±(p 1 ,m)= ^ + i± ^ (2 T | (^1 + m) , v± (p 2)m ) = —-—( J fa- m )\i±) ) (3) 

where f3± = (1 ± /3)/2, and we use the shorthand notation 1^) = \pf). With this definition the massive 
spinors have a smooth massless limit. For m = 0: 

«±(pi) = (1*1 , v±(p 2 ) = |2 T ) . (4) 

We consider now a quark-antiquark vector current. In the helicity basis, we find: 

2J) =M ± ( Pl ,m) 7 ^ T (p 2 ,m) = (^V^) . (5) 

Namely, the vector current has the same functional form as the massless one but replacing the massless 
momenta by the corresponding null vectors used to define the massive spinors. For massive fermions, 
we shall consider also the case of helicity flip along the fermion line, which would vanish in the massless 
limit. We obtain: 
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S»(l±,2f)=u ± ( Pl ,m)Yv ± (p2,m) = -——-(&-&)>> . (6) 

Generalizing these results to the insertion of an arbitrary odd number of gamma matrices between 
the spinors, the following useful properties hold: 

£ W ...m 2 „+i( 1 ± 2 T) = u ± {p 1 ,m)Y 1 ■■■Y 2n+1 v T {p 2 ,m) 

= (3 + {l ± \Y 1 ---l^ n+1 \2 ± )+(3-(l ± \Y 2n+1 ---Y 1 \2 ± ) , (7) 



and 



^i-^n+i^i 2 ±) = u ± ( Pl ,m)7 Ml •■ ■^ 2n+1 v ± (p 2 ,m) 



((2 T | 7 ^ ■ .. T ^n+i|2T) _ (l±| T w .. . 7 w n+ i|!±^ . ( 8 ) 



m 

In the massless limit (3 + — > 1 and /3_ — > 0, and only the first term in Eq. Q survives, while Eq. © 
obviously vanishes. Similar relations can be derived for the insertion of an even number of gamma 
matrices as well, but those are not needed for the current purpose of this paper. 

In the following we will consider the quark-antiquark vector current, with emission of an arbitrary 
number of gluons V 1 — > qq + (n — 2)g. The colour decomposition is such that 

S»(l q ; 2, 3, . . . , n - 1; n q ) = g^ 2 £ (T a *T a * ■ ■ ■ T a " _1 ) S^l q ; 2, 3, . . . , n - 1; n q ) , (9) 

P(2,...,n-X) 

where P(2, . . . , n — 1) is the permutation group of all the gluons, and T ai are colour matrices in the 
fundamental representation. 
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3 Spinorial off-shell currents 



Recursion relations for the calculation of the spinorial current of a quark-antiquark pair and n-gluons, 
where either the quark or the antiquark are off-shell, have been derived long time ago [4J. Those recur- 
sion relations are valid regardless the massive or massless character of the quark-antiquark pair, although 
explicit results where obtained only in the massless approximation. In this Section, we use these recur- 
sion relations to extend their results to the heavy quark case. 

The colour ordered spinorial current of an on-shell quark of four-momentum p 1 and (m — l)-gluons 
of four- momenta p 2 to p m is given in terms of the spinorial current of the on-shell quark with less gluons, 
and the off- shell gluonic current J M of the rest of the gluons: 

m— 1 i 

S(l q , 2,...,m) = -J2 S(l q ; 2, . . . , k)J(k + 1, . . . , m)- , (10) 

where pi >m = p\ + p 2 + . . . + p m and S(l q ) = u(pi, m). For all gluons of positive helicity the gluonic 
current has the form |@]: 

^<><«(« + i)>---0'0 ' 

The null vector £ is the reference gauge vector which is assumed to be the same for all the gluons. 
Conversely, for the off-shell quark current one has 

1 n 

S{m + l,...,n-l;ng) = — J(rn + 1, . . . ,k - l)S(k, . . . ,n - l;n q ) , (12) 

$Wl,n + m k=m+2 

the antiquark carrying four- momentum p n , and S(nq) = v(p n , m). 

It is convenient to work in the gauge where £ = p n , where p n is the reference null vector of the 
antiquark. In that gauge, diagrams with gluons of positive helicity attached to an antiquark of negative 
helicity do not contribute in the massless approximation. This is not the case for massive quarks, but 
that choice of gauge still provides some simplifications. Furthermore: m 2 = For the 

complementary helicity configurations it would be convenient to choose £ = j?i instead. 

The no-gluon, one-gluon, and two-gluon spinorial quark currents are 

-1/2 



S(l+) = u + ( Pl , m) = ^y(m +m) = $T [1| + (£| 



(13) 



Sdt; 2+) = -S{li) it = -S(lt) A mn + m) , (14) 



S(lt-X,S + ) = -(5(lJ;2+)^ + ^(l ? + )7(2+,3+))^— — 

= /co\/a\ — ^ ?i2 3 + — — £ (^i,3 + m) , (15) 

Z/i,3 (£2) (£3) \y 12 (23) / 
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where we defined y l n = p\ n — m . After some algebra we can simplify further the two-gluon current 

S(l» + ;2+,3+) = -|^(JU / dl^W,3 + „,). (16) 
(23)yi2 V(£ 3 ) 2/1,3 / 



Then, we factorize the off- shell current into the form 

S(l 



5'(l+;2+,...,n + ) 



2/l2 



X(l+;2+ 



™ + ) (f|#i,n + m) 



(17) 



where 



^(i 9 + ;2- 



X(l+;2 + ,3- 



J3_ 

(£2) ' 

j_fj2i , a ms\\ 

(23) y li3 ; 



(18) 



The massless limit forces the off-shell massive current to have the following structure 



X(l+;2+,...V 



1 



(23)(34)...((n-l)n) 



(tn) 



+ (3. |£] ,1(1+; 2+ . • • , " H 



(19) 



From Eq. (fT2l . we can easily derive a recursion relation for the X factor of the off-shell current 



*(1^2 + , 



,n+) 



71-2 



#2,n|0 2/12 



(£2)<23)...<(n-l)n) 



+ g 2 + , • • • , fc H , (e(A . + 1)){(A . + + 2)} _ _ _ {{n _ l)n) 



X(l+;2+,...,n-l + )(e|ii w |n] 



(20) 



To simplify this recursion relation further and derive an equivalent recursion relation for the function A, 
we multiply and divide the last term of Eq. d20l) by ((n — l)n), such that 

(f|jfi,»-iM((n - l)n) = (n - l|^i, n -i|0 = (n - l|(^i, n - jrfi,n-2)jfi,r»)IO 

= 1|P?,„ -^l,n-2(^l,n-2 +^n-l,n)|0 

= -Vl,n(£(n - 1)) - (n - %i, n _ 2 +^i, n -2^n-l,ni0 • (21) 

The first term in the rhs of Eq. (f2~Tb generates the first term of the rhs of Eq. ( fl9t . Then, it is clear that 
in order to fulfill the massless limit, the remaining contributions have to either cancel to each other or be 
proportional to /3_. The proof follows by induction. Assuming that Eq. dl~9t is true for n — 1, we show 
that it is also true for n through the recursion relation in Eq. d20l) . For that, we use recursively the Fierz 
identity: 

(k\ = 777T~~ TvT «*(* - l)><ei + (tk)(k - 1|) , (22) 



and 



(£(*-!)> 

(%l,fc +1*1,^+1^10 = (^|2/I,fc-1 + ^l,k-l^k,n\0 ■ 



(23) 
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This results into the following recursion relation for the corresponding A function: 



A(l+;2 + ,...,n + ) ' ^^i+.o+ l-k <*(* + !)> 



Y: A(l+; 2+ . . . , fc + ) fc ^ (e|^ fc+ i,nlO - [2|*,„|0 

(24) 

where 



A(l+;2+) = 0, A(lJ;2+,3 + ) = M. (25) 

2/1,3 



From Eq. (l24h we find 



A(l, + ;2 + ,3 + ,4 + )^(l-2 + ,3 + )(|| + I|^M) . (26) 
To obtain this result we have transformed he last term in the rhs of Eq. (|24h in the following way 

[2|jM6) = A(l+;2+ 3+) [2|yi ^" 10 = A(l+;2+ 3+) + (3|j^ 3 ,nlo) ■ (2V) 

Assuming that for an arbitrary number of gluons larger than two the A function has the form 

2+,. . . , - + ) = 2 + , 3+) (|| + A{lt ' ,2 ^ 3] ' ,n+) ) , (28) 

the A function fulfills the recursion relation 

.4(l + ;2+,..., n +) - 1 r^*„ + ., + ,««# + D> 



Y: A(l+; 2 + , . . . , fc + ) v fc _|_ j (e|^fc + i,nlO " [2|^l«4,n|0 

(29) 



By solving this equation recursively, we get finally 

x, 1+a+ _ ^ 4 A(l+;2+,3 + , 4+) (45)^1^5,^ 10 ~ (£5)[2|^ 23 ^, Wl |Q 
A^ig,/ ,...,n J — 2^ ; 

j = l VlyWl ' ' ' Ul,Wj 

( Wl ( Wl + 1)) • • • K-iK-i + 10 • • • (0^i,^-i^,-i+i,^IO 

(^5)(^i)(0^i + !))••• + 

(30) 

where w k e {5, n}, and w\ < w 2 < ■ ■ ■ < Wj with Wj = n. This completes the calculation of the off-shell 
quark spinorial current. 



3.1 Antiquark off-shell current 

Very similar results can be obtained for the antiquark spinorial current when the antiquark carries nega- 
tive helicity. The no-gluon, one-gluon and two-gluon currents are given by 

s{n- q ) = aT io + p+ 1/2 p| |i] = ~(^n - m)\o , od 
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S(2+;3,) = (^ 3 -m)|£> 



[2| 



{^2)2/23 



S(3 



q J > 



S(2+,3+;4: 



V(£2) 2/2,4 / (23)2/34 



For an arbitrary number of gluons, we have in general 



S(k + , . . . ,n - l + ;n^ ) = (^ fcjn -m)|^)X(A; + ,...,n-l + ;n (7 ) 



2/n— l,n 



with 



X(A; + ,...,n-l + ;nr) 



1 



(*;(* + l))---<(n-2)(n-l)) 
[n - 1| 



where 

A(k + , ...,n-l + ;n~) = A(n -2 + ,n- 1 + ; 



P + [C\A(k + ,...,n-l + ;n^) 



" (e(n-2)) 



+ A(£; + ,...,n- l + ;nr) 



(32) 
(33) 

(34) 



(35) 



(36) 



A recursion relation similar to Eq. (l29l) can also be obtained for this antiquark current, leading to an 
expression similar to Eq. d30b where now the labeling of the four-momenta runs backwards. 



4 On-shell scattering amplitudes 

From the spinorial off- shell current one can calculate the amplitude for the process of n — 2 outgoing 
gluons and a quark-antiquark pair. This amplitude is obtained from the quark-gluon spinorial current 
with n — 2 gluons, by removing the propagator of the off-shell quark, contracting the current with the 
antiquark spinor and imposing momentum conservation 



5(1+; 2+, . . . , n - 1 + ; nf) = -^1 <?(1 + )X(1+; 2 + , . . . , n - 1 + ) (£\S(nJ) 
h i yi2 y 



Pl,n=0 



The helicity conserving amplitude annihilates, and for the helicity flip amplitude we get 

m 3 2/lin _ 1 A(l+;2 + ,...,n-l+) 



S(l + ;2+,...,n-l+;n+) 



9 ' ""' ' 9y /3+yia (23><34>...((n-2)(n-l)>(nl> 

In particular 



2/l,ra-l=0 



(37) 



(38) 



s « ;2+ ' 3+;4 ^^(Hiy <39) 
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5 Quark-antiquark vector current 



From the results obtained in the previous sections we can also calculate the off-shell quark-antiquark 
vector current 

n-l 

2,...,n-l;n g -) = £ S(l q ; 2, . . . , k)-fS(k + 1, . . . , n - 1; n 9 ) . (40) 



For all gluons of positive helicity all the terms contributing to the sum in Eq. (l40l) have the form 

2 + , . . . , n - 1+; nj) cc + m) 7 ^ fc+ i,n - m))^) . 

Then, after some algebra we find that 

S£(l+;2+ ...,n-l + ;n-) oc -m(£\jk ,«-i7l£) , 



, ...,n- 1 ,nT 



oc (e|j<i,fc7^*+i,n|e]-2m 2 ^ 



(41) 

(42) 
(43) 



Particularly interesting is the helicity conserving current. In that case all the terms are proportional to 
the same quantity, where now the quark four-momentum has been replaced by its reference null vector 
pi. Since in Eq. (l42l only massless vectors appear some properties that hold for massless scattering 
amplitudes can be generalized to the massive case directly without further proof. For some applications 
see Section|6j We obtain the following expression for the helicity conserving vector current: 

1 / /3+[12] /Ml(n-l)] 



^(l+;2+,...,n-l + ;nr 



+ 



m 



(£1) 



2/12 



,2 n-2 



(23><34>.-.<(n-2)(n-l)> \{t{n-l))y 12 
1 



,n- 1 + ) + 



2/n— l,n 



-A{V 



— £r £<*(* + !)> 

2/12 2/n-l,n [If J k=2 



[12] 



+ /5_[ie]A(l+;2^ 



X 



[(n-l)l] 



+ /3+[l£)A(k + r 



n 



'(£|#l 1 n-lTlO- 



1)) ' - 

In the massless limit, the current in Eq. (l44h agrees with the Berends-Giele current 0|: 

(f|*i,»-i7"IO 



5W;2 + ,...,n-l + ;nr) 



In the following we present some simple examples. 



5.1 Single gluon emission 



=0 (12)(23)---((n-l)0 



(44) 



(45) 



Let us consider a quark-antiquark vector current with emission of one single gluon off the quark 
antiquark pair. For the helicity conserving vector current, we find: 

f/3+ 0\ [12] 



V2/12 2/23/ (£2) 
The helicity flip configuration is also quite simple, we get 

S"(l+;2 + ;3}-) = ro 
which clearly annihilates in the massless limit. 



12710= /3^%|#i27l0 



2/12 2/23 



'^-^m2<ft-ft)" + JUi(*- + %i2> 

,2/12 2/23/ (£2) (fl) V2/12 2/23/ 



(46) 



(47) 
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5.2 Two gluons vector current 



The helicity conserving vector current with emission of two gluons is 

1 f [12] j3 + /3_[13]\ 



S"(l+;2+ 3 + ;4, 



m 



[21] [31] 



(23) V (£3) Vu 2/34 m J 2/12 2/34 (£2) (£3) [If] 



! [32] 



1 



+ 



(23) (fl) I 2/12 1/1,3 2/34 2/2,4 



(48) 



6 Higgs production together with heavy quarks 

One of the main objectives of the LHC program is the discovery of the Higgs boson. Gluons couple to 
Higgs through a top quark loop dominantly, where the top can be integrating out leading to an effective 
interaction. MHV rules for Higgs plus multiparton amplitudes have been presented in Refs. lETl H7\ 
in that effective theory. The naive off-shell continuation of the Higgs vertices depends of the gauge 
spinor however, and therefore do not work as MHV amplitudes. This problem is solved by splitting the 
Lagrangean into a self-dual and an anti-self-dual part, and introducing a complex scalar where the Higgs 
is one of its components = h(H + iA). 

In the heavy top quark effective theory amplitudes of 0-scalars coupled to a quark-antiquark pair and 
any number of gluons of positive helicity vanish ll22l 

M&Qiigt, •••,<?;UO = o. (49) 

In this section we shall extend this result to the heavy quark case. The proof follows almost straitforward 
from the amplitudes that we have calculated in Section The helicity conserving vector current of a 
heavy quark-antiquark pair and m — 2 gluons is 

£"(1+ 2+ . . . , m - 1 + , n- q ) = iV 2 , m _! (£|#i,m-i7l£> , (50) 

and the vector current of all gluons of positive helicity |4| reads 

J"(m+ . . . ,n - 1+) = A^-i^lTVr^-ilO . (51) 

The normalization factors N are irrelevant for the present proof. What is important is that Eq. d50l 
has the same functional form in the massive and in the massless case, where p x is the null vector used 
to define the quark spinor. The 3-point <pgg vertex of two off-shell gluons of momenta k x and k 2 and 
Lorentz indices \i\ and /x 2 is given by ll22l : 

^/w/ia^ 1 ' ^2) — 2//J1/J2 ki ■ k 2 — ^i M 2^2/ji + *£/ii/uai/ii/2&i 1 &2 ! ' 

Then, one can compute A n {<P, Qi, gt \ • • • > flvt-i) Qn ) ^y contracting the vertex with the currents in Eq. d50l) 
and Eq. (ED, where k x = p 1>m -i + p n = pi, m -i + p n , and k 2 = p m , n -i- 

The 3-point <pgg vertex contracted with the quark-antiquark and the Berends-Giele currents gives 

V*,S»r = (h ■ k 2 )(S ■ J) - (ki • J)(k 2 ■ S) 

1 71—1 

+ ((^T^iii] + s»j» ■ 

i=m 
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Fierzing the gluonic current 



ra-1 



J=2N m>n - 1 £ + (jO, (54) 



]=m 



and after some algebra, we find 



V+S"J" = N m ,^ 1 %- (55) 



where ^i] is the commutator of the quark-antiquark current with its total momentum. Fierzing the 
quark-antiquark current 



m— 1 



?=2AWE(I*K£I + I0W)<6'>, (56) 



the 3-point vertex obviously cancels since in Eq. (l55l) : = and = 0. 



7 Conclusions 



We have calculated several tree-level multigluonic helicity amplitudes with heavy quarks. Our approach 
is based on the Berends-Giele recursive relations. We have shown that the method can be easily applied 
for the evaluation of multipartonic scattering amplitudes with heavy particles, and leads to compact 
results. We have considered scattering amplitudes with all gluons of the same helicity only, but other 
helicity configurations can be calculated efficiently as well. Furthermore, interesting phenomenological 
implications can be derived from the general structure of multipartonic amplitudes. 
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A Spinors and polarization vectors 

Massive spinors of four-momentum p M and mass m can be constructed with the help of two null vectors 
p M and as follows: 

u±(p,m) = —^L—{^ + m) 1^) , u± (p, m) = ^ (q T | (J + m) , 

v±(p,m) = (fi T \^±\ ($ ~ m ) 1^) » v±(p,m) = j^r^{q T \ (tf - m) , (57) 
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where N = (p • q/p ■ q)^ 1 ^ 2 is the normalization factor. Those spinors trivially fulfill the Dirac equation: 

{jf> — m) u(p, m) = , {jf> + m) v(p, m) = , (58) 
and the completeness relation 

m -a^a = i> + m , XI v -^\ =i>-m. (59) 
The polarization vector of an outgoing gluon of momentum k^ is defined as 

^■« = ± Sp- (60) 

where £ M is an arbitrary gauge reference momentum satisfying £ 2 = 0. The polarization vector contracted 
with the gamma matrices reads as follows: 

Other practical properties helping to simplify calculations are 

(^1^,0^=0, 1jf t (k,Z)\i k ) = 0, (62) 

and 

(ri^±(^O = 0, ^±(A;,Oie ± ) = 0. (63) 
Finally, if we choose the reference momenta for all gluons to be the same, then it follows that 

e ± (k l ,O-e ± (k r O = 0. (64) 
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